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Abstract 


I  have  developed  a  graphical  user  interface  (GUI)  that  makes  the 
process  of  evaluating  the  radar  range  equation  faster  and  more 
convenient.  From  this  GUI,  one  can  enter  the  input  parameters 
necessary  to  compute  the  signal-to-noise  ratio  (SNR)  and  probability 
of  detection  (Prf)  of  a  radar  system  as  a  function  of  range,  and  view  the 
results  graphically.  The  GUI  and  the  SNR  and  Prf  computations  are 
coded  using  the  MATLAB  programming  language.  I  review  the 
derivation  of  the  equation  used  to  evaluate  the  P(/  of  a  signal  detected 
by  a  linear  quadrature  detector  and  provide  examples  of  the  GUI's 
operation. 
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1.  Introduction 


Evaluating  the  radar  range  equation  and  obtaining  probability  of  detec¬ 
tion  (Pd)  as  a  function  of  range  are  frequently  useful  when  performing  a 
radar  systems  analysis.  A  MATLAB  graphical  user  interface  (GUI)  has 
been  developed  to  do  this  conveniently  and  flexibly.  Shown  in  figure  1, 
the  GUI  allows  parameters  to  be  entered  from  either  the  pop-up  boxes, 
sliders,  or  editable  text  boxes.  The  values  taken  from  the  text  boxes  are 
used  to  calculate  curves  of  signal-to-noise  ratio  (SNR)  versus  range  from 
the  radar  range  equation  (fig.  1  (a)),  and  these  curves  are  plotted  as  a 
function  of  range  for  three  values  of  the  radar  cross  section  (RCS).  A 
second  plot  of  Pa  (fig.  1  (b))  is  generated  based  on  the  computed  SNR 
values  and  a  probability  of  false  alarm  (Pfa)  entry.  This  evaluation  tool, 
consisting  of  a  MATLAB  GUI  and  MATLAB  code  for  computation,  is 
intended  to  be  used  as  a  first  cut  for  system  design,  since  it  does  not  have 
the  level  of  detail  to  support  an  advanced  study.  Although  it  is  flexible 
enough  to  support  most  surveillance  radar  designs,  such  as  the  Army 
Research  Laboratory's  (ARL's)  low-cost  enabling  radar  technology 
(LCERT)  radar  or  the  Yuma  Proving  Ground  (YPG)  instrumentation  radar 
programs,  it  would  have  to  be  modified  to  address  other  type's  of  radars. 
However,  it  would  not  be  difficult  to  extend  or  modify  the  capabilities  of 
this  tool  or  to  add  features,  such  as  a  detailed  system-loss  budget  includ¬ 
ing  atmospheric  attenuation  or  antenna  shapes  other  than  circular.  Modi¬ 
fications  could  also  be  made  to  this  MATLAB  code  to  change  the  param¬ 
eter  computations  or  to  read  out  intermediate  results. 
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2.  Graphical  User  Interface 

The  GUI  is  called  with  the  command  "rrefig"  at  the  MATLAB  prompt. 

The  command  "help  rrefig"  will  print  a  synopsis  of  how  to  use  the  GUI. 
To  use  the  GUI,  enter  the  desired  parameter  values  into  either  the  sliders, 
pop-up  boxes,  or  editable  text  boxes.  After  each  entry,  a  new  set  of  SNR 
and  Pd  curves  is  computed.  Since  this  computation  takes  about  3  s  (almost 
entirely  because  of  the  Pd  computation),  the  response  of  the  sliders  may 
be  somewhat  slow.  The  values  entered  into  the  pop-ups  or  sliders  are 
copied  to  the  associated  editable  text  boxes.  Conversely,  a  value  entered 
into  an  editable  text  box  causes  the  associated  pop-up  to  read  "see  edit 
box."  However,  a  value  entered  into  an  editable  text  box  does  not  cause 
any  change  in  the  position  of  the  associated  slider.  The  transmitted  band¬ 
width  and  window  function  inputs  are  not  used  as  inputs  to  the  radar 
range  equation.  Rather,  they  are  used  to  compute  the  theoretical  range 
resolution  that  would  be  obtained.  The  example  computation  shown  in 
figure  1  is  based  on  a  conceptual  design  for  a  projectile  tracking  instru¬ 
mentation  radar.  The  RCS  values  used  are  based  on  a  range  of  generic 
artillery  projectiles.  Another  example,  shown  in  figure  2,  is  based  on  the 
LCERT  system  design.* 


*Eric  Adler  et  al,  "Low-Cost  Enabling  Technology  for  Multimode  Radar  Requirements,"  Record  of  the  IEEE 
1998  Radar  Conference,  12-13  May  1998. 
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Figure  2.  LCERT  field 
test  parameters: 

(a)  signal  to  noise 
versus  range  and 

(b)  probability  of 
detection  versus 
range. 
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3.  Radar  Range  Equation 


One  can  use  a  form  of  the  radar  range  equation  found  in  Skolnik's  hand¬ 
book*  to  obtain  SNR  as  a  function  of  range  and  RCS: 


SNR  = 


PtGA1aX2GihwGl 


rdop 


(4  7rykTBR4LN 


f 


(1) 


where 


Pf  =  peak  transmitted  power, 

Ga  =  antenna  gain, 

<7  =  radar  cross  section  of  target, 

A  =  wavelength  of  transmitted  carrier, 

G,to,  =  time-bandwidth  product  gain, 

Gdop  =  Doppler  processing  gain, 
k  =  Boltzman's  constant, 

T  =  ambient  temperature, 

B  =  modulation  bandwidth  of  transmitted  signal, 

R  =  range  to  target, 

L  =  system  losses,  and 
Nf  =  noise  figure. 

Although  the  equation  is  exact,  the  parameters  used  to  obtain  the  SNR  are 
often  approximations  or  based  on  engineering  judgment.  Obviously,  the 
user  must  apply  knowledge  of  the  problem  and  experience  to  obtain  a 
reasonable  result,  and  the  value  of  this  tool  is  that  it  makes  this  process 
more  efficient.  The  details  of  the  SNR  computation  may  be  found  in  the 
comments  in  the  MATLAB  code  shown  in  the  appendix  of  this  report.  The 
antenna  gain  is  computed  based  on  entering  the  diameter  of  a  circular 
antenna  with  an  efficiency  of  60  percent. 


*M.  Skolnik,  Radar  Handbook,  2d  edition,  1990,  p  2.6. 
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4.  Probability  of  Detection 

In  this  section,  I  will  discuss  the  derivation  of  the  equation  that  I  used  to 
compute  the  Pd  as  a  function  of  SNR  for  a  given  Pfa,  The  derivations  are 
based  on  course  notes  by  Trunk*  and  a  book  by  McDonough  and 
Whalen.+ 1  will  consider  the  situation  where  an  input  signal  v(t)  is  proc¬ 
essed  through  a  linear  quadrature  detector,  as  shown  in  figure  3.  The 
integral  function  in  the  quadrature  detector  represents  a  low  pass  filter. 
The  output  of  the  quadrature  detector  is 


q  = 


v(t)sin(coct)dt\  + |  J  v(t)cos(coct)dt 


1/2 


(2) 


where  r=  measurement  period,  0<t<z,  and  (0C  =  carrier  frequency. 

One  can  now  postulate  two  possibilities  for  v(t).  It  may  consist  only  of 
white,  Gaussian,  zero-mean  noise,  or  it  may  consist  of  a  sinusoidal  signal 
plus  the  noise.  Thus  I  define  two  hypotheses: 

H0:  v(f)  =  n(t) 

Hp  v(f)  =  n(t)  +  A  sin( coct  +  6)  ,  (3) 


where  n(t ) 


coc 

A 

e 


G( ) 


G( 0, 

carrier  frequency, 
amplitude  (peak-to-peak), 
random  receive  phase, 
variance  of  noise  signal, 
Gaussian  (normal)  distribution. 


and  the  signal-to-noise  ratio  of  power  is  defined  as 


SNR  = 


(4) 


Figure  3.  Linear 
detector. 


*G.  V.  Trunk,  Detection  and  Estimation  Theory  Handbook ,  notes  for  Johns  Hopkins  course  525.728,  p  60. 

+r  N.  McDonough  and  A.  D.  Whalen,  Detection  of  Signals  in  Noise ,  2d  edition.  Academic  Press,  1995, 
pp  159-166. 
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(In  eq  (4),  I  define  the  SNR  as  power  rather  than  as  energy  as  in 
McDonough's  definition  (eq  (7.32)  and  (7.54)).  If  ris  normalized  to  equal 
1  in  equation  (7.32),  then  the  power  for  a  complex  signal  is 


P=  2 


iJl 


+{i%) 


A ‘ 
2 


(5) 


where  the  factor  of  2  outside  the  brackets  is  introduced  for  convenience, 
as  it  was  in  equation  (7.25)  of  McDonough's  work. 

The  Neyman-Pearson  criterion*  is  typically  used  in  radar  detection 
problems  to  determine  an  expression  for  Pd  for  a  given  Pfa.  It  specifies  a 
procedure  to  test  one  hypothesis  against  another  and  is  a  most  powerful 
test — Pd  is  maximized  for  a  Pfi  less  than  or  equal  to  some  value.  To  sim¬ 
plify  the  derivation,  one  can  first  assume  continuous  probability  density 
functions  (pdf's)  for  the  distributions  of  interest.  Then  one  can  assume 
that  the  envelope  of  the  input  signal  does  not  vary  during  the  measure¬ 
ment  period  t.  If  it  had,  one  could  still  show  that  the  test  is  uniformly 
most  powerful  for  the  amplitude  and  could  proceed  the  same  way.  One 
could  then  compute  the  Pd  for  a  target  with  Rayleigh  fading,  or  for  a 
Swerling  II  target  when  multiple  coherent  processing  intervals  (cpi's)  are 
integrated  noncoherently.  Finally,  one  assumes  that  although  the  phase 
angle  6  of  the  input  signal  is  unknown,  it  has  a  uniform  random  pdf. 
Although  the  random  variable  in  the  signal  makes  H1  a  composite  hy¬ 
pothesis,  it  can  be  integrated  out,  with  no  change  in  the  form  of  the  pdf 
for  Hlr  and  now  Hj  becomes  a  simple  hypothesis. 

One  then  defines  the  Pfa  as  the  probability  that  a  detection  Dj  is  declared 
when  hypthesis  H0  in  equation  (3)  is  correct: 


Pfa  =  piD  1 1  Ho)  =  £  Po(q)d(l  /  (6) 

where  T  is  defined  as  the  detection  threshold  level  and  po(q)  is  the  pdf 
associated  with  hypothesis  H0.  Similarly,  the  Pd  is  defined  as  the  probabil¬ 
ity  that  a  detection  Di  is  declared  when  the  hypothesis  Hi  in  equation  (3) 
is  correct: 

Pd  =  P{D1 1  Hj)  =  £  p-fq)dq  ,  (7) 

where  p^q)  is  the  pdf  associated  with  hypothesis  Hl.  The  threshold  T  will 
be  some  function  of  q  and  is  found  by  using  the  likelihood  ratio  test.  This 
test  states  that  you  decide  Hi  if 


Hq)  = 


pM 

Po(<?) 


>T  . 


(8) 


*R.  N.  McDonough  and  A.  D.  Whalen,  Detection  of  Signals  in  Noise,  2d  edition,  Academic  Press,  1995 
pp  159-166. 
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To  find  the  ratio  in  equation  (8),  one  must  determine  the  numerator  and 
denominator  pdf  s.  McDonough*  derives  the  pdf  of  a  linearly  detected 
narrowband  signal  plus  narrowband  noise.  It  is  known  as  the  Rician 
density  function,  and  for  this  problem  is  given  by 


where  I0  =  Bessel  function  of  the  first  kind  order  zero.  If  the  signal  q  is 
equal  to  zero,  then  I0(0)  =  1,  and  the  resulting  density  will  be  that  of 
integrated  noise,  which  is  the  Rayleigh  probability  density  given  by 


(10) 


The  threshold  can  now  be  found  by  substituting  equations  (9)  and  (10) 
into  equation  (8): 


In  equation  (11),  only  the  term  inside  the  argument  of  the  Bessel  function 
contains  the  term  q.  So  the  threshold  test  reduces  to 


I 


^  T  , 


(12) 


where  constants  have  been  absorbed  into  qt •  Since  the  Bessel  function  is  a 
monotonically  increasing  function  of  q,  the  ratio  test  can  be  performed  as 
a  linear  function  of  q  so  that 

q  >  T  .  (13) 

From  equation  13,  one  can  set  T  =  qt  and  use  equations  (6)  and  (10)  to  find 
qt  in  terms  of  a  given  value  of  P/«: 


\dq  =  exp 


(14) 


By  rearranging,  one  obtains 

^  =  [-21n(f»]'/2-  <15> 


*r.  N.  McDonough  and  A.  D.  Whalen,  Detection  of  Signals  in  Noise,  2d  edition.  Academic  Press,  1995, 
pp  129-130. 
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Substituting  equation  (9)  into  equation  (7),  the  Ptf  is 


H,(#xp 

The  substitution,  v  =  q/c rfi,  is  now  made  in  equation  (16)  to  put  the 
integral  in  the  form  of  the  Marcum  Q  function,*  which  is  more  convenient 
for  calculating.  Also,  by  using  the  expression  for  SNR  in  equation  (4), 
equation  (16)  becomes 


V  +  A2) 


207 


Aq 


°l  a2 


(16) 


q(\/ snr,^) 


-q,/a 


'exp 


■  ( v2  +  SNR) 


I0(sfSNRv)dv  .  (17) 


This  is  the  expression  used  to  compute  the  Pd  for  a  given  Pfa,  where  qt/o 
is  found  from  equation  (15).  The  code  used  to  evaluate  equation  (17)  in 
terms  of  equation  (15)  is  found  in  the  appendix.  It  is  important  to  note 
that  the  SNR  as  defined  in  equation  (4)  is  lower  than  the  definition  in 
Skolnik+  by  a  factor  of  2  or  3  dB.  The  Skolnik  definition  is  used  in  the 
MATLAB  GUI  so  that  the  SNR  calculated  from  the  radar  range  equation 
matches  his  curves  in  figure  2.3  of  his  publication  rather  than 
McDonough's  curves  in  figure  7.3  of  his  publication.  Also  note  that  the 
SNR  as  defined  above  is  in  units  of  volts  per  volt,  not  decibels.  In  equa¬ 
tion  (1),  the  computation  is  performed  in  decibels.  So  one  converts  to  the 
form  used  in  equation  (17): 

'/SNR  =  l°l°gio(^fp)  ■  (18) 


*R.  N.  McDonough  and  A.  D.  Whalen,  Detection  of  Signals  in  Noise ,  2d  edition,  Academic  Press,  1995, 
p  131. 

+M.  Skolnik,  Radar  Handbook ,  2d  edition,  1990,  p  2.6. 
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5.  Conclusion 


A  MATLAB  GUI  has  been  developed  that  can  be  used  to  quickly  and 
easily  determine  the  SNR  and  Pd  as  a  function  of  the  radar  range  equation 
parameters  and  the  Pfa .  I  have  used  it  to  plot  curves  of  SNR  and  Pd  versus 
range  for  two  different  radar  examples:  the  YPG  instrumentation  radar 
and  the  LCERT.  The  equation  used  to  compute  Pd  as  a  function  of  SNR 
and  Pfa  was  derived  using  the  Neyman-Pearson  criterion  and  procedure. 
One  could  further  develop  the  capabilities  of  this  GUI  by  adding  new 
functionality  and  features. 
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Appendix.  Code  for  Radar  Range  Equation  and  Pd  GUI 


%  FILE:  rre 

%  DATE:  7/19/98 

% 

%  This  function  is  the  target  of  callbacks  in  the  rrefig  GUI.  Radar  range  equation 
%  parameters  and  the  Pfa  are  read  from  the  editable  text  windows,  and  used  to 
%  compute  SNR  and  Pd  as  a  function  of  SNR  for  the  given  Pfa.  The  vulues  from  the 
%  popups  and  sliders  are  read  into  the  associated  editable  text  boxes.  After 
%  reading  the  parameters,  the  SNR  as  a  function  of  range  is  computed  for  3  RCS 
%  values.  Then,  the  Pd  as  a  function  of  range  is  computed  based  on  the  SNR  values 
%  from  the  radar  range  equation.  See  technical  report  "A  MALAB  Radar  Range  Equation 
%  and  Probability  of  Detection  Evaluation  Tool" . 

%  The  MATLAB  copy  figure  function  does  not  copy  the  popups,  sliders,  or  boxes. 

%  To  copy  everything  in  the  GUI  to  a  eps  file  with  a  tiff  preview  feature,  use 
%  command  "print  -depsc2  -tiff  <path . f ilename>" . 


%  trans_pwr 
%  antenna_diam 
%  rf__loss 
%  frequency 
%  rcs_start 
%  noise_figure 
%  pulsewidth 
%  cpi_pulses 
%  window 


%  Transmitter  power  in  Watts. 

%  Antenna  diameter  in  meters 
%  RE  system  hardware  loss 

%  Transmitter  frequency  in  GHz  (1  GHz  =  1) 

%  RCS  in  dBsm 

%  Receiver  LNA  noise  figure  in  dB 
%  Pulsewidth  in  microseconds. 

%  Number  pulses  in  coherent  processing  interval 
%  Rectangular  for  no  window,  or  Hamming 


%  snr 


%  Power  signal- to-noise  ratio,  in  dB. 


%  action  ==  1  reserved  for  initial  conditions  if  required, 
function  rre (action) 


%************  READ  PARAMETER  VALUES  FROM  GUI  ***************************** 


%  Execute  this  block  from  every  action  on  the  GUI  callbacks.  For  popups 
%  and  sliders,  their  callbacks  put  their  values  into  the  editable  text  block 
%  which  is  then  read  in  this  section, 
if  action  ~=  1 

h  =  f indob j  (  'Tag'  ,  'Freqe');  %  Finds  the  object  with  the  tag  'Freqe'  . 

frequency  =  eval(get(h,  'String7));  %In  this  case,  eval  is  like  str2num. 

%  frequency  =  eval (get (gcbo ,  'String'))  used  in  Celsius  example 

%  Use  code  below  to  get  value  from  popup;  see  p  3-22  of  GUI  manual 
%h  =  f indob j { ' Tag ' ,  'Freq'); 

%value  =  get(h,  'Value'); 

%string  =  get(h,  'String'); 

%frequency  =  str2num ( string(value) ) ; 

%  If  the  action  is  from  the  freqe  editable  text  block,  set  the  value  of  the 
%  freq  popup  to  "see  edit  box".  Similarly  for  the  other  popups. 
if  action  ==  3 

h  =  f indob j ( 'Tag ' ,  'Freq'); 

set(h,  'Value',  5)  %  Set  the  value  of  the  popup,  'Freq',  to  5,  which  is 

%  "see  edit  box"  in  the  string  cell  for  'Freq' 

end 


u 


Appendix 

h  =  f indob j ( 'Tag' ,  'Power e'); 
trans_pwr  =  eval (get (h,  'String')); 
if  action  ==  4 

h  =  f indob j ( ' Tag' ,  'Power'); 
set (h,  'Value' ,  6) 

end 

h  =  f indob j ( 'Tag' ,  'Diame'); 
antenna_diam  =  eval (get  (h,  'String' ) ) ; 

h  =  f indob j { 'Tag' ,  'RCSe'); 
rcs_start  =  eval (get (h,  'String')); 

h  =  f indob j ( 'Tag' ,  'Losse'); 
rf__loss  =  eval  (get  (h,  'String')); 

h  =  f indobj ( 'Tag' ,  'Nfe'); 

noise_figure  =  eval (get (h,  'String')); 

h  =  f indobj ( 'Tag' ,  'Pwe'); 

pulsewidth  =  eval (get (h,  'String')); 
if  action  ==  5 

h  =  f indobj ( 'Tag' ,  'Pw'); 

set (h,  'Value ' ,  7 ) 

end 

h  =  f indobj ( 'Tag' ,  'Numbere' ) ; 
cpi_pulses  =  eval (get (h,  'String')); 
if  action  ==  6 

h  =  f indobj (' Tag' ,  'Number'); 
set (h,  'Value' ,  7) 

end 

h  =  f indobj ( 'Tag' ,  'Pfae'); 
pfa  =  eval (get (h,  'String')); 
if  action  ==  7 

h  =  f indobj ( 'Tag' ,  'Pfa'); 

set (h,  'Value' ,  8) 

end 

h  =  f indobj ( 'Tag' /  ' Strt_rnge ' ) ; 

start_range  =  eval (get (h,  'String'))  *  1 . e3 

h  =  f indobj ( 'Tag' ,  'Stp_rnge'); 
stop_range  =  eval (get (h,  'String'))  *  1 . e3 ; 

h  =  f indobj ( 'Tag' ,  'BWe'); 

bandwidth  =  eval(get(h,  'String')); 
if  action  ==  8 

h  =  f indobj ( 'Tag' ,  'BW' ) ; 
set (h,  'Value ' ,  6) 

end 

h  =  f indobj ( 'Tag' ,  'Window'); 
value  =  get (h,  'Value'); 

string  =  get(h,  'String'); 
window  =  string {value} ; 


end 
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EVALUATE  RADAR  RANGE  EQUATION  ******************* 


c  =  3e8 ;  %  Velocity  of  light  (m/sec) . 

A  =  pi  *  antenna_diam  ^2  /  4; 
wavelength  =  c  /  ( frequency  *  1 . e9 ) ; 

antenna_gain  =  10  *  loglO  (.6  *  4  *  pi  *  A  /  wavelength  A  2); 


four_pi  =  10  *  loglO ( ( 4  *  pi)A3); 


pt  =  10  *  loglO  (trans_pwr)  ; 

lambda_sq  =  2  *  10  *  loglO (c  /  (frequency  *  l.e9)); 
ktb  =  10  *  loglO (1.38e-23  *  300  * (bandwidth  *  l.e6)); 
t_bw__gain  =  10  *  loglO  (pulsewidth  *  bandwidth); 
dop„gain  =  10  *  loglO <cpi_pulses ) ; 

range  =  1  inspace  (star t__range,  stop_range,  100); 
rcs(l)  =  rcs_start; 
delta_rcs  =  -10; 

for  ii  =  1  :  3  %  Loop  over  3  RCS  values 

snr(ii,  :)  =  pt  +  lambda_sq  +  2  *antenna_gain  +  t_bw_gain  +  dop_gain  +  rcs(ii)... 
-  four__pi  -  ktb  -  40  *  loglO  (range)  -  noise_figure  -  rf_loss; 
rcs(ii+l)  -  rcs(ii)  +  delta_rcs; 

end 


%  *****************  EVALUATE  RANGE  RESOLUTION 


*************************** 


res  =  c  /(2  *  (bandwidth  *  l.e6)); 
if  window  ==  'hamming 
res  =  res  *  1.44; 

end 

h  =  f indob j ( 'Tag' ,  'Res'); 

set(h,  'String',  res) 


%  *****************  EVALUATE  Pd  FUNCTION  ******************************** 

%  Taken  from  McDonough  and  Whalen,  "Detection  of  Signals  in  Noise",  2nd  ed 
%  eqn  7.53 

global  alpha  %  Make  alpha  global  to  pass  it  to  marcum_q_fn. 

beta  =  sqrt ( -2* ( log (pf a) ) )  ; 

%snrl  =5  :  15; 

rangel  =  1 inspace (star t_range,  stop_range,  50); 
lrl  =  length(rangel); 
for  ii  =  1  :  3 
for  jj  =  1  :  lrl 

%  In  the  Marcum  Q  function,  the  SNR  is  defined  as  a  voltage  ratio,  but 
%  the  rre  SNR  is  a  power  ratio.  To  take  the  square  root,  divide  by  20  instead  of 

10. 

alpha  =  10 . A ( (snr ( ii ,  2  *  jj)  +  3)  /  20); 

%  Add  3dB  to  SNR  to  match  Skolnik,  p2.20 
%  Don't  add  3dB  to  SNR  to  match  McDonough,  p  2  61. 

%  Care  should  be  taken  in  the  use  of  tolerance  for  the  integration.  Larger  toler¬ 
ances 

%  will  run  faster,  but  will  generate  errors  in  Pd  for  small  values  of  Pd.  The 
errors 

%  are  not  noticeable  in  the  linear  scale  plots. 
pd(ii,  jj)  =  1  -  quad  ( 'mar  cum_q_fn' ,  0,  beta,  [l.e-1  l.e-3]); 

end 
end 
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11  =  sprintf { 'RCS  =  %3.0f  dBsm' ,  rcs(l))?  %  Use  these  in  the  legends. 

12  =  sprintf ('RCS  =  %3.0f  dBsm',  rcs(2)); 

13  =  sprintf ('RCS  =  %3.0f  dBsm',  rcs(3)); 

haxes  =  f indob j ( 'Tag ' ,  'Axesl'); 

set (haxes ,  'nextplot',  ' replacechildren ' ) ;  %  These  two  lines  keep  plots  from  printing 

%  to  outer  figure  box. 

axes (haxes ) 

x_lims  =  get  (haxes,  'xlim')  ;  %  Use  this  to  make  xaxis  in  second  plot  consistent 

plot (range,  snr(l,  :),  'k-\  range,  snr(2,  :),  range,  snr(3,  :),  'k-') 

grid  on;  title (' Signal-to-Noise  vs.  Range') 

xlabel ( 'Range  (km)');  ylabel (' Signal  to  Noise  Ratio  (dB)  ' ) 
legend (11,  12,  13,  0) 

haxes  =  f indob j ( 'Tag' ,  'Axes2'); 

set (haxes,  'nextplot',  'replacechildren'); 

axes (haxes) 

plot (rangel ,  pd(l,  :),  'k-',  rangel,  pd(2,  :),  'k-.',  rangel,  pd(3,  :),  'k-') 
axis  (  [x_lims  ( 1)  x_lims(2)  0  1])%  %  Use  x  axis  limits  from  first  plot, 

grid  on;  title  (' Probability  of  Detection  (Pd)  vs.  Range  (Skolnik  p2.20)') 
xlabel (' Range  (km)');  ylabel ('Pd') 
legend (11,  12,  13) 


%  FILE:  marcum_q_f  n .  m 

%  DATE:  7/29/98 

%  Use  this  function  as  integrand  to  evaluate  the  Marcum  Q 
%  function  that  provides  the  probability  of  detection  for 
%  a  single  pulse  out  of  a  quadrature  detector.  Called  from 
%  rre.m  which  calculates  the  radar  range  equation. 


function  f  =  marcum_q__fn  (v) 
global  alpha 

f  =  v  .*  exp ( - (v. A2  +  alpha. ^2) /2) .*  besseli(0,  alpha. *v) ; 
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